Space-Time Description of Scalar Particle Creation 
by a Homogeneous Isotropic Gravitational Field 
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Abstract. We give the generalization of the method of the space-time description of particle creation 
by a gravitational field for a scalar field with nonconformal coupling to the curvature. The space-time 
correlation function is obtained for a created pair of the quasi-particles, corresponding to a diagonal form 
of the instantaneous Hamiltonian. The case of an adiabatic change of the metric of homogeneous isotropic 
space is analyzed. We show that the created pairs of quasi-particles in de Sitter space should be interpreted 
as pairs of virtual particles. 
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1. INTRODUCTION 

Quantum field theory in curved space-time is 
nowadays a sufficiently deeply elaborated area of 
theoretical physics (see the monographs [HEI) with 
important applications in cosmology and astro- 
physics. In particular, creation of particles with 
GUT-scale masses by the gravitational field of the 
early Universe may be used for an explanation of 
the observed visible and dark matter density |3j. 

In the description of particle creation by the gravi- 
tational field, a widely used method is that of instan- 
taneous Hamiltonian diagonalization [T] suggested 
by A. A. Grib and S. G. Mamayev [4]. A detailed 
study of the created pair formation process was per- 
formed in [5] for the case of a scalar field conformally 
coupled to the curvature. The method of a space- 
time correlation function suggested in [5] made it 
possible to distinguish real created particles from 
virtual ones, to reveal the role of horizons in par- 
ticle creation etc. 

In quantum field theory in curved space-time, one 
frequently considers scalar field nonconformally cou- 
pled to gravity, in particular, the minimally coupled 
one. In such cases, the contributions related to non- 
conformal coupling may be dominant both in the 
particle creation effect [6] and in the vacuum aver- 
ages values of the stress-energy tensor (see, e.g., [7]). 
Such a coupling with the curvature also takes place 
in the case of massive vector mesons (the longitudi- 
nal components [I] ) . Conformal invariance is lacking 



in the graviton equations as well Therefore, it 
appears to be necessary to generalize the method 
of studying the particle creation process suggested 
in ^ to the nonconformal case. Moreover, some 
authors [3] consider the absence of such a general- 
ization as an argument in favor of choosing only the 
conformal coupling in the wave equation for a scalar 
field. 

The present paper suggests a generalization of the 
space-time description of particle creation 5] to the 
case of nonconformal coupling. In [Section "2} we 
perform quantization of a general-type nonconfor- 
mally coupled scalar field in homogeneous isotropic 
space. In [Section 3[ we carry out diagonalization of 
the generalized Hamiltonian built in [lOj . which, in 
the nonconformal case, allows one to solve the well- 
known problem of an infinite density of quasiparti- 
cles created [11 . In lScction 41 we build a space-time 
correlation function for quasiparticles corresponding 
to the diagonal form of the instantaneous Hamilto- 
nian and study the case of an adiabatically changing 
metric. In IScction"5l we consider particle creation 
and the space-time correlation function in de Sitter 
space. The IconclusionI brieflv sums up the results of 
the paper. 

We use the system of units in which h — c = 1 . The 
signs of the curvature tensor and the Ricci tensor are 
chosen in such a way that Rik = R^. 
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where r*^, are the Christoffel symbols. 
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2. SCALAR FIELD IN CURVED 
SPACE 

Consider a complex scalar field f{x) with mass m, 
the Lagrangian 



L{^) = Wd,^*dkV - (m^ + Vg) , (1) 
and the corresponding equation of motion 



(V'V, + Vg + m^) ^{x) = , 



(2) 



where Vi are covariant derivatives in TV-dimensional 
space-time with the metric gik, g — det{gik), Vg is 
a function of invariant combinations of the metric 
tensor gn^ and its partial derivatives. 

Eq. ([2|) is conformally invariant if m = and 
Mj = ^cR, where R is the scalar curvature and 
^c = iN-2)/[A{N- 1)] (conformal coupling). The 
case Vg = corresponds to minimal coupling. An 
arbitrary Vg leads to the advent of third- and higher- 
order derivatives of the metric in the metric stress- 
energy tensor of the scalar field and consequently in 
the Einstein equations. 

It is well known that additional terms with higher- 
order derivatives appearing in equations lead to rad- 
ical changes in the theory even if the coefficients of 
these terms are small. 

If one requires that the metric stress-energy ten- 
sor should not contain derivatives of the metric of 
orders higher than two, then the following function 
is admissible as Vg: 



where 



R. 



dct 



GB 



Vg=£,R + CRlB, 
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R' 



(3) 



(4) 



(the Gauss-Bonnet coupling [T2]). 

Let us note that for TV = 4, with constant (p{x), 
the contribution to the metric stress-energy tensor 
from the term with Rq^ is absent because the cor- 
responding variation derivative vanishes !13j. But 
for a variable ip{x), a contribution from such terms 
could be taken into account if the constant ( with 
the dimension (mass)~^ is nonzero. 

Accounting for a possible coupling between a 
scalar field and the Gauss-Bonnet invariant Rq^ 
may play an important role in the early Universe; 
effects from a nonzero value of the parameter C in 
scalar field equations may appear in black-hole ra- 
diation, may affect the parameters of the so-called 
boson stars etc. The question of the values of the 
parameters C and ^ are ultimately related to the area 
of the experiment. 



Furthermore, without specifying Vg, let us con- 
sider an iV-dimensional homogeneous isotropic 
space-time, choosing the metric in the form 

ds^ = g.kdx'dx'' = a^{T]) {di]^ - df) , (5) 

where dP — ja/idx^dx^ is the metric of an {N — 1)- 
dimensional space of constant curvature K = 0, ±1. 

The complete set of solutions to Eq. ^ in the 
metric ([5]) may be found in the form 



(p{x) 



where 



ip{x) 



a(7V-2)/2(^) 



g'i{v) + ^'{v)9xiv)^0, 
n^r^)^{m^ + Vg-^,R)a^ + \^ 



(6) 

(7) 

(8) 



N-2 



Uj(x), (9) 



the prime denotes a derivative with respect to 
the conformal time 77, and J is the set of indices 
(quantum numbers) numbering the eigenfunctions 
of the Laplace-Beltrami operator Aat-i in (A^— 1)- 
dimensional space. 

According to the Hamiltonian diagonalization 
method U (see the case of an arbitrary function Vg 
in [S]), the functions g\{ri) should obey the follow- 
ing initial conditions: 

g'xivo) = i ^ivo) g\{vo) , \9x{vo)\ = fi^'^^'^ivo) ■ 

(10) 

To perform quantization, let us expand the field 
(p{x) in the complete set of solutions (O 



(p{x) = / d^l{J) 



(+) A+) 



(-) A-) 



(11) 



where d^{J) is a measure on the set of quantum 
numbers, 

(12) 

and require that the standard commutation relations 
hold for a"^^ and a^^\ 

Let us build the Hamiltonian as the canonical 
one for the variables (^(x) and if>*{x), for which the 
equation of motion does not contain their first-order 
derivatives with respect to the time r] [14]. Recall 
that the equations of motion do not change after 
adding a full divergence d,P jdx^ to the Lagrangian 
density L(x). 
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Let us choose, m the coordinate system (fy, x), the 
vector 

(J*) = (/7c(^*^(iV-2)/2,0,...,0), 

where 7 = det(7c,^), c — a' /a. Then, using the 
Lagrangian density L'^{x) = L{x) + dJ^/dx^, we 
obtain for the momenta canonicaUy conjugate to (p 
and (p*: 



connected with a 7^', a"^' by time-dependent Bo- 
goHubov transformations: 



dL 



A 



dif' 



dL 



A 



d(p*' 



Vl^'. (13) 



respectively. Integrating the Hamiltonian density 
h{x) = lp'tt + (^*'7r* — L^{x) over the hypersurface 
S: 7y = const, we obtain the foUowing expression for 
the canonical Hamiltonian: 



(14) 



(m2+ Vg) a' - (2c' + {N- 2^) 



(see a justification of such a choice of the Hamilto- 
nian in [inilll] and in Section 3). 

The Hamiltonian (|14p may be written in terms of 

the operators Qj and Uj in the following way: 



dn{J) 



E 



(-)„(-) 



(15) 



where 

E.i = 



\9'x? 



(16) 

and we have chosen such eigenfunctions $,/(x) that, 
for arbitrary J, there is such J that <i>}(x) = 
i?j$j(x), = 1 , (J = J, = i?j). Such a choice 
is possible due to completeness and orthonormality 
of the set <i>j(x). 

In spherical coordinates of a homogeneous 
isotropic space, if J = {A, to}, we have J — 
{A, -to}, dj = {-ir (see [I]). 



3. HAMILTONIAN 
DIAGONALIZATION 

The Hamiltonian (fT5|) will be diagonal at the time 
instant r/o with respect to the operators aj , a j , 
under the conditions pn]) . Diagonalization of the 
Hamiltonian at an arbitrary time instant 77 is carried 

out in terms of the operators h j (jf) and hj '{r]), 



=«}('/) &J ^(r?)~/5,/(ry)^j6;f^(r]). 



*(-) 



(+V, 



(17) 



where the functions aj{rj) = aj{r]) and f3j{rj) = 
Pj{ri) satisfy the initial conditions |aj(77o)| = 1, 
/3j(?7o) = and the identity \aj{Tj)\^ - |/3j(r?)|2 = 
1. Substituting p?| and the conjugate expres- 
sions to (I15p . one can obtain an expression for the 
Hamiltonian having the same form (|15p but with the 



replacement d 



(±) „(±) 



¥f\bf^ and 



^ bEj = Ej{\aj\^+ 2Re(Fjaj/3}7?}), 

(18) 

i^j ^ t^^j = -2aj/3ji?j£;j + a2^^j +/3}7?jFJ. (19) 
From the diagonality requirement for the Hamil- 
tonian with respect to the operators bj (jf) and 
b'-j'\r-i) at the time 77, i.e., bFj{rj) = 0, (for n^{7]) > 
0), it follows: 



.9a 



where Xj = Xj is an arbitrary complex function of 
time with a unit absolute value. Therefore further 
we will use the operators 



(21) 

which, due to P?|) and (^0]) . do not depend on the 
specific choice of the functions Xj iv) ■ The operators 
Cj {71) and Cj (r/) obey the same commutation re- 
lations as and a[^\ 

Substitution of ([2T|) . taking into account (fT7|l and 
their conjugate expressions, into ([TS]) leads to the 
following expression for the Hamiltonian: 



(-)/ 



i/(77)= dfi{J)n{r^)(^c 



(+)„(-) 



+ 0^4+)). (22) 



Thus the energy of quasiparticles corresponding to 
the diagonal form of the Hamiltonian ([TS]) is equal 
to the oscillator frequency i^{ri) (unlike the Hamil- 
tonian built from the metric stress-energy tensor of 
a nonconformal scalar field f6j). 

Using (III), mi), jini), and one can verify that 
the operators c[^\'r]) obey Heisenberg's equations of 
motion: 



dr/ 



^' ,q(Tl)„(T) 



(23) 
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The second term in the right-hand side of is 
connected with re-definition of the particle notion 
at each time instant. 

An expansion of the field operator (p{x) in the 

dfi{J) 



operators Cj^ 



/2n 



<f,Kx)4+)(r,) + <i>,Kx)4 ^(,7) 

(24) 

follows from lU), (HH), and (HOD- The equa- 

tions of motion hold for each mode in ((24)) separately. 

Consider the question of particle creation in a non- 
stationary metric. We suppose that the quantized 
scalar field is in the state |0), annihilated by the op- 
in the vacuum state for the 



(—) * ( 
erators a\ \ 



I.e. 



V ' "'J 

instant 779. At the time rj, the vacuum state is the 

state 



|0^), defined by the equalities 



„(-) 



(25) 



The state |0) contains, at the time instant 77, 
|/3j(?7)P pairs of particles and antiparticles corre- 
sponding to the operators c[^\'r]) (see [T|). The den- 
sity of the created particle pairs may be calculated 
(for the quasi-Euclidean metric with K = 0) by the 
formula 



n{ri) = 



B 



N 



(26) 



where Bn = [2^-3^(^-i)/2r((iV- l)/2)] \ T{z) is 
the gamma function, S\(ri) = |/3a('7)P (in a homoge- 
neous isotropic space, = |/?a|). For = 4 and 
K = 0,-1, for the number density of the particle 
pairs created, the following formula is valid (see [1]): 



71(77) 



1 



27r2a3 



Sxiv) A2 dX, 



(27) 



For K = 1 (spherical space), the set of eigenfunc- 
tions of the Laplace-Beltrami operator A^r-i is dis- 
crete, and the formula for the number density of the 
created particle pairs A^ = 4 has the form (see [1]) 



1 



27r2a3 ^ 

A=l 



5a(77)A2 



(28) 



Using ([20]) and that the function 

is a first integral of the equation ([7]), equal to — 2i 
according to the initial conditions (jlOp. we obtain: 



Sxiri)^^{\9'x? + ^'\9>. 



(29) 



As it was shown in [TU], Sx{ri) ~ A^^ as A — > 
00. Therefore, in four-dimensional space-time, the 
number density of particles created, defined by the 
Hamiltonian (jl4p diagonalization method, is finite in 
the nonconformal case as well. Let us note that a di- 
vergent expression for the number density of created 
nonconformal scalar particles, obtained with another 
choice of the Hamiltonian in Ref. [TT] , has been one 
of the reasons for a criticism of the Hamiltonian di- 
agonalization method as a whole in the well-known 
book 121. 



4. THE SPACE-TIME 
CORRELATION FUNCTION 

To study the space-time characteristics of the cre- 
ated quasiparticles, we apply the approach suggested 
in [S]. We use the notion of a particle's localized 
state introduced by Newton and Wigner [T5]. By 
analogy with the Newton- Wigner operator for a free 
field (see, e.g., [TB]), we introduce creation operators 
of a localized state of a particle and an antiparticle: 



'^(+)(77,x) =a- — (77) /d^(J)<i>}(x) 0^^(77), 



Af-l 



:(+)/ 



^(+^(77,x) =a-— (77) yd/i(J)$}(x)c7^(7;).(30) 

By analogy with the case of a conformal scalar 
field considered in |1|, we introduce the operator of 
particle number in the volume V 

Nv^J d^-'xV^^ ^<+^(77,x)^^r^(77,x), 

V 

(31) 

where (^"i^g = det ((^"^^ga/s), and 

is the induce metric tensor on the hypersurface rj = 
const. Using (fTT)) . ((2T|) . pO|) and the properties of 
the eigenfunctions fE'j(x) (see, e.g., 1 , §9.1), it can 
be shown that the expression for the number density 
of the particles created 



N~l 



( + )/ 



71= {Q\Nv\Q)/V 



(32) 



obtained with the aid of ((3T|) . reproduces Eqs. ([261) - 
(Ell). 

As in [5] , we consider as a characteristic of the spa- 
tial distribution of the quasiparticle pairs created, 
the matrix element 



i?o(?/,x,x') 



(0,1^^1 ^(77,x) ^[ \v,^')\0) 
(0,|0) 



, (33) 
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which has the meaning of the probabiHty amplitude 
that a quasiparticle created is located at the point x 
at the time instant rj while the antiquasiparticle is 
at the point x'. Using ^ and §0^, we 

obtain: 



i?o(??,x,x') = - 



1 



iv) 



rf/i(J)$j(x)$}(x')PAW, 

(34) 



where 



Pxiv) - mgx ' g;^)/mgx + g'x) ■ (35) 



The function P\{ri) satisfies the following equation 
and initial condition: 



Px{m)^0- (36) 



Consider the case that the metric is changing adi- 
abatically: 







n' 




< 






fi2 



< 1. 



(37) 



Furthermore, we denote M = - >?la. In the 
general case, M — M(r\). For a conformally cou- 
pled scalar field, M = m. U M — const, the con- 
ditions §7i) hold if h(t)/M'^ < h/M < 1, where 
h{t) = d{t)/a is the Hubble parameter. 

If the conditions ([57]) and a'(?7o) = are satisfied, 
an approximate solution to (|36p has the form 



(38) 



Let us find an expression for the space-time cor- 
relation function (j33p in the approximation con- 
sidered in the quasi-Euclidean metric (i.e., ds^ = 
a?{ri){drf' — dx^dx"")). The eigenfunctions of the 
Laplace operator Ajv_i in the coordinates x" are 

$j(x) = (2^)(i-^)/' eM-^^ax") , (39) 

Jla^l = Conse- 



where — cx3 < Aa < 
quently. 



J {A— const) 



(40) 



where I] = A <5(| A | - A), p |x - x'|, 
Ji/(z) are Bessel functions. Substituting ((38)) . 
into (|34[) and performing integration, we obtain 

i{Ma)^' { M 



and 



i?0(»/,X,x'): 



^(^)^ X^(Mr), (41) 



where r — pa and i^i/(z) are MacDonald's functions. 



If M — const, which is the case for a conformally 
coupled scalar field and for an arbitrary coupled field 
in de Sitter space, from ((4T|) we obtain 



i?o(i,x,x') 



-iKPhjt) / M \- 
87r2 V2^/ 



KN^{Mr). (42) 



Let us further consider spatial sections with X = 
±1. For A'^ = 4, the space-time metric may be writ- 
ten in the form 



df = dx^ + fix) {M^ + sin^^dif^) 



(43) 



where /(x) = sinh(x), X, sin(x) for K ~ — 1,0,+1, 
respectively. Meanwhile, 

V $}(x)<i>,;(x') = A^^. (44) 



J (A— const) 



2^2 /(^) ' 



where p is the geodesic distance between the points 
X and x' (see [1]). Therefore, from ((34| and (|38p we 
obtain in the hyperbolic {K = —1) case: 



i?o(»y,x,x') 



-z(Ma)2' p 
IGTT^a^ sinhp 



i^o(A^r-). (45) 



If |z| > 1, then K^{z) ~ 1/^(2^) e-^ Conse- 
quently, if the metric is changing adiabatically, the 
function |i?o| in the quasi-Euclidean (see (|4T|) ') and 
hyperbolic space-time decreases exponentially with 
growing p at spacings r 3> 1/Af, i.e., exceeding the 
Compton wavelength. 

In spherical case, the spacing p changes in the 
range < p < tt. A substitution of ([38]) and p4)) 
into (IMll leads to 



i?o(ry,x,x') 



-i(Afa) 



2/ 



A sin Ap 



167r2a3 sinp (Af2a2 + ;^2)3/2 



(46) 

The sum in the right-hand side of (|46|) , after a chain 
of transformations, may be presented in the form 



"S^(p+27rn) A'o (Ma(p-f 27rn)) 



n=0 



n=l 



(27rn-p) Kq (Afa(27rn-p)) . 



(47) 



From the asymptotic properties of the function 
Kq{z) it follows that, for p < 1 and Ma > 1 (the 
distance between particles od a pair and a particle's 
Compton wavelength are much smaller than the cur- 
vature radius of space), in this representation one 
could retain only the term p KQ{Map), and there- 
fore Eq. ((46|) takes the form 



i?o(»y,x,x') 



~i{Ma 



i2/ 



IQir^a^ sin p 



Ko{Mr) . (48) 
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Therefore, in the spherical case |i?ol also decreases 
exponentially with growing p at distances r ^ 1/Af, 
exceeding the Compton wavelength. 

Thus if the metric of a homogeneous isotropic 
state is changing adiabatically, the space-time cor- 
relation function _Ro(?7,x, x') is exponentially small 
for r ^ \/M. This indicates that the corresponding 
quasiparticles are virtual pairs with the character- 
istic correlation length equal to \/M. Real particle 
creation is exponentially small and does not mani- 
fest itself in perturbation theory. 



n/M^ 



0.0015 
0.001 
0.0005 



1 



M/H 



5. PARTICLE CREATION 
IN DE SITTER SPACE 

Let us consider de Sitter space and take the metric 
in the form (HI) with K = Q and 



a = aie 



Ht 



1 



(49) 



t e (-00, -f oo) r] e (-0O, 0). 

Solutions to Eq. ([7]) with Vg = and the condi- 
tions ([To]) for rjo ^ —oo have the form 

9xiv) = eii-''i/(2)(-A^)e"" , (50) 

(2) 

where Hi (z) is a Hankel function, ao is an arbi- 
trary real constant, 

^ = ^^-(§)', A/ = + (51) 

and R = N{N- l)H^ . 

Furthermore, assuming m? + (^ — ^c)R > 0, from 
(I26l). (I29ll and (EOl we obtain 



N 



(52) 



i.e., the created particle number density in de Sitter 
space is time-independent! The result of a numerical 
computation for the function Fn{M/H) in the case 
= 4 is represented in Fig. [T] 

At M/H ^ 1, the metric is changing adiabatically 
and, as shown in Section 4, real particle creation 
does not occur. 

In the general case, substituting the exact solu- 
tion ([501) into ([Ml) and ([35]) and using (gO]), we ob- 
tain 



M^-i ■' 1' H 



(53) 



Figure 1: The number density of quasiparticles cre- 
ated in de Sitter space. 

|Ro|/M^ 
0.01 




r-M 



2 3 4 5 6 7 

Figure 2: The correlation function for M/H = 0.01. 



where r = a(i)|x' — x|, i.e., the correlation func- 
tion, expressed in terms of r (the "physical" dis- 
tance between the quasiparticles in a pair) is time- 
independent. 

Examples of numerical calculations for A^ = 4, 
M/H = 0.01 and M/H = 0.2 are given in Figs. [2] 
and [21 respectively. 

In both cases, the correlation function decreases 
exponentially at distances between the quasiparti- 
cles exceeding the Compton wavelength Ic = 1/M. 

Thus real particles creation in de Sitter space does 
not occur. The quasiparticle pairs being created, 
whose density has been calculated and shown in 
Fig. [U should be interpreted as pairs of virtual par- 
ticles. 

As has been noticed in [T] , the absence of real par- 
ticle creation in de Sitter space is confirmed by the 
local nature of the vacuum stress-energy tensor and 
by vanishing of the imaginary part of the effective 
Lagrangian. 
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|Ro|/M^ 




— r-M 

2 4 6 8 

Figure 3: The correlation function for M/H = 0.2. 

6. CONCLUSIONS 

In this paper, for a scalar field nonconformally 
coupled to the curvature, we gave a generalization 
of the method of space-time description of particle 
creation by the gravitational field. In a homoge- 
neous isotropic space, we have introduced the cre- 
ation operators ([50)1 of localized one-particle states 
and the operator (|3ip of particle number in a speci- 
fied volume. We have obtained the expressions 
and ([55]) for the space-time correlation function 
of a pair of created quasiparticles corresponding to 
a diagonal form of the instantaneous Hamiltonian. 
We have analyzed the case of adiabatic changes in 
the metric of a homogeneous isotropic space. The 
expressions (gl]), (gH), (gn]) and (gH]) have been ob- 
tained for the correlation function of a pair of quasi- 
particles created. It has been shown that the corre- 
lation function exponentially decreases at spaces ex- 
ceeding the Compton wavelength, and consequently 
real particle creation is suppressed. Particle creation 
in de Sitter space has been considered, and, from, 
the properties of the space-time correlation function 
for a pair of quasiparticles created, it has been con- 
cluded that such a pair should be interpreted as a 
pair of virtual particles. 
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